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Abstract 

Classical solutions of the Yang-Mills-dilaton theory in Euclidean space-time are investigated. 
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of solutions exist in finite region of parameter space and discuss this issue in detail in different 
dilaton field normalization. 
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I. INTRODUCTION 



Bartnik and McKinnon (BK) discovery of static, spherically symmetric, globally regular 
solutions of the Einstein- Yang-Mills (EYM) theory [l| has opened new avenue of research 
in General Relativity. Soon after BK discovery related (colored) black holes were found and 
different generalisations were investigated (see e.g. Q for review and references). Interesting 
modification of pure EYM theory is inclusion of a dilaton field, which is predicted from many 
modern theoretical considerations like string theory and extra dimensions. So, EYM-dilaton 
theory was investigated [1, 0, 0] and solitons and black holes similar to BK case were found. 
It came as a surprise that even eventually simpler YM-dilaton system in flat space-time has 
infinite tower of globally regular solutions similar to the BK case. These solutions can 

be labelled by an integer n - number of nodes of gauge field amplitude. The lowest solutions 
in EYM and (E) YM-dilaton theory have one node and the whole tower of solutions turns 
out to be unstable with n th solution having 2n unstable modes: n in 'gravitational' sector 
and n in 'sphaleron' sector 0, 10, [ll|, 12 



Note that most of these investigations were performed in four dimensional space-time with 
Lorentzian signature. In such case it was found that electric part of non-abelian 577(2) gauge 



field should necessarily vanish for asymptotically flat solutions [13j, [14(. Situation is changed 



in Euclidean sector: electric field here plays role similar to the Higgs field in Lorentzian 



sector. Euclidean version of the EYMD theory was investigated [15| and solutions with 
non- vanishing electric field and non- vanishing non-Abelian charges were found recently. 

The aim of present paper is to study Euclidean solutions in the YM-dilaton theory. The 
rest of the paper is organized as follows: in next Section we formulate our model, derive 
equations of motion and asymptotic behaviour of finite action solutions. In Sect. Ill we 
discuss some special solutions of our system. In Sect. IV we present our numerical results. 
Sect. V contains concluding remarks. 



II. THE MODEL 



Starting point of our investigation is the Euclidean version of Yang-Mills theory coupled 
to a dilaton field <b: 



St 



la,. 



l <j>+ 1 ^F^jd 4 x, (1) 

where F^ u is the non-Abelian gauge field strength and 7 and g denote the dilatonic and 
gauge coupling constants, respectively. 



A. The Ansatz and equations 



In this study we will restrict ourselves with the SU(2) gauge g roup and will be inter- 
ested in spherically symmetric solutions. Following Witten [ly, [l7| the general spherically 
symmetric SU(2) Yang-Mills field can be parameterized as follows: 



Al = X -u{r) , 



A,- (-jak-Ek 



1 - W(r) 



+ 



XjX a 



A\ XjX a 
— + -A^A 2 



(2) 



with r = y xf being the radial coordinate, i, j = 1, 2, 3 are spatial indices and a,b — 1, 2, 3 
group indices. 



2 



The reduced action S re d is denned as follows: 



Tma 



S E = dTS re d , (3) 







with t = xo being Euclidean time. Without loss of generality we can put A\ = A 2 = 0. The 
reduced action then reads: 



'red 



p 27</> (1 _ U/2\2 i 



where the prime denotes the derivative with respect to r. 

The equations of motion which follow from the reduced action (jlj) read: 

W= -2^'W'- (l -^ 2)W + Wn\ (5) 

7/' 2 

u" = -2- - 2 7 0V + 4^ 2 « , (6) 



The rescaling 



2 v . 2 ^\^2 (l-^ 2 ) 2 rV 2 , 

2r 2 

— > — , r — > — r , it — > — m , (8) 

7 3 7 

removes the gauge and dilatonic coupling constants g and 7, respectively, from the equations 
of motion. We thus set without loosing generality 7 = g — 1 in the following. 



B. Asymptotic behaviour 

Close to r = we find a 3-parameter family of solutions regular at the origin, which can 
be parameterized in terms of b, u\ and <po: 

W(r) = l-br 2 + 0(r 4 ) , (9) 
u(r) = u x r - yr 2 + 0(r 3 ) , (10) 

<f>{r) = 0o + 2e 2 * (b 2 + ^f)r 2 + 0(r 3 ) . (11) 

The behaviour of the solutions at infinity is more involved. Assuming a power law be- 
haviour for u and 0, for leading behaviour of W we obtain expressions containing the Whit- 
taker functions, which in leading order decay exponentially. We thus obtain: 

W{r) = Ce- u °° r r Q °(l + 0(-)) , (12) 
u(r) = Uoo -^ + 0(\), (13) 

4>(r) = 0OO- — + 0(i). (14) 
So, at infinity there are five free parameters 0oo, Uoq, Q e , Qd and C. 
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Note that the shift of dilaton field 0^0 + for any finite value can be compensated 
by a rescaling of r and u as follows: 

r^re^, u -> we"^ . (15) 

This symmetry can be used to normalize the solutions e.g. as 0(r = 0) = or 0(r = oo) = 0, 
and reduces the number of free parameters in the system by one. 



III. SPECIAL SOLUTIONS 
A. Analytic solutions 

The equations of motion (jHJ), (jSJ), (|7|) have three important analytic solutions. The first 
one is the vacuum solution with 

W{r) = ±1 , u{r) = , 0(r) = . (16) 

The second one is a simple, dyonic-type generalization of dilatonic monopole solutions 
which have been found previously in the case for u(r) =00]: 

W{r) = , u(r)=u , 0(r) = ln(^j— ) , (17) 

where u is an arbitrary constant. As in the pure magnetic case this solution can be obtained 
by integrating the first order Bogomol'ny-type equation as follows. The reduced action (jl]) 
for W(r) = and u(r) = const can be written as 

Sred = \ J dr (V - + e*|§° • (18) 

Thus for configurations with values 0(0) = — oo and <f>(oo) = the action ([TBI has a minimum 
S re d — 1 and attains this minimum for solutions of the first order Bogomol'ny-type equation 

r 2 0'-e </, = O. (19) 

Solutions of this equation automatically satisfy the equations of motion ©-([7]) with W = 
and u = const. Integrating (jT9l) we obtain the logarithmic behavior ffTTl) for the dilaton 
field. 

The third solution is well known Bogomol'ny-Prasad-Sommerfield (BPS) monopole solu- 
tion, which in our notations reads: 

W = — T — - , u = coth(r) - - , = 0. (20) 
smh(r) r 



B. Solutions for u = 



The system of Eqs. (EE) for u = was first analyzed in 0] and it was shown that a series 
of solutions exist indexed by the number of nodes n (with n > 1) of the function W. In the 
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normalisation 0(0) = each of these solutions can be uniquely characterized by a value of 
the 'shooting' parameter b entering in (Q; for the first few n we have 

61 « 0.26083 , b 2 « 0.35352 , 63 « 0.37500 . (21) 



The precise numerical values of parameters 6 n for n = 1, ..,8 can be found in [8j. The 
solutions with a higher number of nodes are very well approximated on an increasing domain 
of r by 

W(r) = -^sin ( ^-^lnr + 5 J , u(r) — , 0(r) = lnr, (22) 



which is solution of the linearized equation of W 



W" + -W + \w = , (23) 



with W\ and 5 arbitrary constants. The limiting solution with infinitely many zeros of W 
described by the asymptotic behavior (I2"2"j) has b — > ps 0.37949. 

The values of masses (reduced action) of static solutions in the normalisation </>(oo) = 
lie between Mi ~ 0.80381 for the solution with one node and = 1 for the limiting 
solution 0, |s| . 



C. Related Bessel function 

In the case of a non-zero electric field, the nature of limiting solution is changed. For 

u = q and = lnr, (24) 
for the gauge amplitude W in the linear regime, W « 1, we get equation 

W" + -W + (\ - q 2 )W = , (25) 
whose solution is a spherical Bessel function 

W(r) = W hessel = ^K^(qr) , (26) 
V r 2 

where W2 is an arbitrary constant. 

The function Wbessei has the following properties, where for convenience we use x = qr. 
For i>1 the solution decays according to W = exp(— x)(l/x + 0(l/x 2 )). Close to the origin 
the solution is not analytic. Solving the equation numerically, we see that the for x < 1 
the solution develops oscillations whose amplitudes become larger while the limit x — >• is 
approached. This is illustrated on Fig. [TJ The first nodes appear for 

X! w 0.0367 , x 2 w 0.001 , x 3 w 0.000026. (27) 

These zeros are separated by extrema, alternatively maxima and minima. The first of these 
points appear at 

xi « 0.1265 , x 2 ^ 0.0032 , x 3 w 0.0008. (28) 
where, up to an arbitrary normalisation, W takes values 

W(zi) = 1 , W(£ 2 ) ~ -6.12 , lV(5 a ) w 35.56 (29) 



5 




1e+1 



FIG. 1: The solution of the Bessel equation (|25p . W and W' are plotted as functions of x. 



IV. NUMERICAL RESULTS 



A. Solutions in the gauge 0(oo) = 

We have solved the above equations numerically using the differential equations solver 
COLSYS 18] . To construct these solutions, we use the following boundary conditions: 

W(Q) = 1 , u(0) = , 0(r)'| r =o = (30) 



at the origin and 

W(oo) = 
at infinity. The parameters b res , u" 



tt(oo) = Mqo , 0(oo) = (31) 

u r ^ s and mass M res for the solutions in the normal- 
ization = are related to the parameters b, %, and mass M in the normalisation 
0o — with a simple shift and rescaling: 







res 




J<f>o 



b , u 



res 
1 



'u 



OO ; 



M T 



(32) 



Practically, u(oo) is the parameter that is fixed by hand in the boundary conditions. 
From a numerical solution with fixed n and tt(oo) one can extract the values of U\ and b, 
and other parameters which enter in asymptotic behaviour Eqs. (I MTTT) and Eqs. fll2H14p . 

We studied in details the solutions of the above equations with zero, one and two nodes 
of the gauge field function W(r). Examples of such solutions are given in Fig. [20 for 
u res (oo) = 0.1 respectively for the gauge field function W(r) for zero, one and two nodes, 
and the corresponding profiles for the functions u and 4>. 



B. Pattern of the solutions 



In order to relate the pattern of solutions available for u ^ to the solutions of we 
convert our solutions to the normalization <po = by fl32l) . The values of the two "shooting 
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r 



FIG. 2: The gauge field function W is shown for zero, one and two node-solutions, respectively, 
for u res (oo) = 0.1. 




FIG. 3: Same as Fig|2]for the gauge field function u r 



parameters" b and u\ then fully characterize the solutions and describe a point P(b, Ui) on 
a plot of initial data. On the U\ = 0-axis, the numerical values found for b agree with the 
solutions of 0] and leads to the points (b n , 0) where b n are given with a high accuracy in 

0-0- 

The relations between parameters b , Uoo and ui > are shown in Figs. EH Let us first 
discuss FigIS] where the values of b are plotted as functions of the parameter for the first 
few values of n. We see immediately that a branch with no node of the function W exist. In 
the limit — > the function W converge to the constant value W = 1 on this branch and 
we recover vacuum solution Eq. (fl6l) . When is increasing, b clearly tends to a constant 
value for 14(00) — > u(oo) cr . The actual value of this constant depends on the number of 



nodes of the gauge field function and increases for increasing number of nodes [19j . Let us 
point out that when Uoo becomes large, the solution appears very concentrate in the region 
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FIG. 4: Same as Fig(2]for the dilaton field function ft 



of the origin and the numerical analysis becomes very difficult. 

Complementary to Fig. \5\ in Fig. [6] we show the dependence of the parameter u\ on 
the parameter u{oo). The parameter u\ first increases with u(oo), reaches a maximal value 
and then decreases and reaches a constant value for -u(oo) — > u(oo) cr . Again for a given 
value of u(oo), the value of u\ increases with increasing number of nodes of the gauge field 
function. Finally, in Fig. [7] we summarize the pattern of the '"initial data"' of our solutions 
by presenting the dependence of the parameter 6 on u\. The three curves finish at 

(b, Mi) ^(0.3687,0.2763) for n = , 
(6,ui) ^(0.3796,0.00707) for n = 1 , 

(6,«0 ^(0.3776,0.00056) for n = 2 . (33) 

On this figure we also see in particular the following properties: (i) For < b < b\ there 
exist only solutions with no nodes of W having "monopole asymptotics" , W^oo = 0. 
An example of such a solution is given by (b, u\) = (0.2,0.2693). (ii) For values of b, 
b\ < b < b 2 there exist solutions with zero and one node. Fixing for example b = 0.27, they 
are characterized by the two points = (0.27,0.30222) and (6,1*0 = (0.27,0.00237) 

and correspond to solutions with monopole asymptotics having zero nodes of W and one 
node of W respectively, (iii) For values of b such that b 2 < b < b% there exist solutions with 
zero, one and two nodes. Taking e.g. b = 0.36, the three points (6, u\) = (0.36,0.28366), 
(6, Mi) = (0.36,0.015866) and (6, Ui) = (0.36,0.000247) with the same parameter b define 
solutions with monopole asymptotics having respectively zero, one and two nodes of W. 

In order to illustrate the critical phenomenon which appears in the large Uoo limit we 
refer to FigJHland Fig. [9] where the profiles of the one- node solutions are presented for three 
values of u^. Several properties of the solutions can be observed on these figure. Namely for 
increasing we observe that : (i) the function u(r) increases more and more strongly in 
the region of the origin and reaches its asymptotic value when W is still very close to W = 1; 
(ii) the dilaton function becomes very close to ft(x) = ln(r) in an increasing interval of r 
where W(r) << 1. The figures further suggests that the function W(x) approaches closer 
and closer the profile of the Bessel function between the second extremum and infinity. 



0.4 
0.35 
0.3 
0.25 
.Q 0.2 
0.15 
0.1 
0.05 




10 



10" 
U(°°) 



zero nodes 

- one node 

- two nodes 



FIG. 5: The dependence of the parameter b on u(oo) is shown for zero, one and two node solutions. 
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FIG. 6: Same as Fig|5]for the parameter u%. 

The evidence of this statement is made particularly apparent on Fig JH] where solutions are 
superposed for three different values of Uoo and where the corresponding function Missel is 
represented by the bullets line. We see in particular that the profiles of W (represented for 
— 0.5, 1, 10 respectively by the dotted, dashed and solid lines) become closer and closer 
to the bullets curve representing the Bessel function (the curve on Fig. 1 has of course to 
be appropriately renormalized and the argument x has to be rescaled). This rather good fit 
between the function W(r) presenting n-nodes and the function Wbessei breaks down in the 
region of the origin (the function Wb eS sei{x) indeed becomes very singular for x — > 0). 

To compute the reduced action of the solutions, it is convenient to return to (fi(oc) = 
gauge. For small values of the following behaviour of the reduced action A = S re d is 
found 



A « 8.9^ for n = , A w 0.803 + 0.5^ for n = 1 , A w 0.965 + 0.1 Uoo for n = 2 . (34) 
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FIG. 7: The plot of initial data is given for the zero, one and two node solutions. 

More details of the evaluation of the action for zero-node solutions are shown on Fig. [TD1 
Within our numerical accuracy, we find that A — > 1 in the critical limit as its expected from 
the property of limiting solution Eq. (|T7|) . 



1.4 




1e-5 1e-4 1e-3 1e-2 1e-1 1e+0 1e+1 1e+2 

r 

FIG. 8: The profiles of the gauge function W and u for one node solutions are presented for 
Uqo = 0.5, 1.0, 10.0 respectively by the dotted, dashed and solid lines. The Bessel function is 
represented by the bullets. 
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FIG. 9: The profiles of the gauge function W and the dilaton for one node solutions are presented 
for Uqo = 0.5, 1.0, 10.0 respectively by the dotted, dashed and solid lines. 

C. Solutions in the gauge u(oo) = C 

When the equations are analyzed in the 'gauge' such that <f)(oo) = and increasing the 
parameter Uoo by hand, the solutions become rapidly very concentrated in a region around 
the origin and the numerics becomes hard, as mentioned already above. As an alternative to 
understand the critical phenomenon limiting the solutions for a maximal (but n-depending) 
value of u(oo), we studied solution with fixed li^ and increasing by hand the value 0(0) as a 
boundary condition. The value 0(oo) can then be determined numerically. Different values 
for can be connected by means of 

u^Au , r^j , 0^0- log(A) . (35) 

We will discuss in details the results in this gauge for the solutions with zero and one nodes of 
W. We would like to point out that no new solutions were discover in the present approach, 
it is rather a complimentary way of describing the solutions obtained with the normalisation 
0(oo) = 0. 
Zero node for W. 

Let us first discuss the solution corresponding to zero node of the function W. If we solve the 
equations with a fixed value for and changing the value of 0(0) by hand, the numerical 
results show that when 0(0) decreases, the values 0(0) and 0(oo) become very close to 
each other and that the function <j>(x) becomes practically constant. As a consequence the 
equations for W(r) and u(r) reduce to the equations of the monopole of SU(2) gauge theory 
in the BPS limit. This correspond to the first part of zero node branch labelled '1' on FigfTOl 
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The fact that the parameter b and u\ tends to zero in the limit 0(0) — > — oo in the same 
limit on Figs. 013 is due to the fact that the natural radial variable of the monopole solution 
is rescaled by a factor exp(— 0(0)), so that the monopole solution appears 'diluted' by the 
dilaton in the corresponding limit. For small enough 0(0) this is confirmed by the numerical 
results which show clearly that the profiles of W, u hardly differ from the (suitably rescaled) 
BPS solution Eq. (1201. 

When we increase the value 0(0), it appears that solutions can be constructed up to a 
maximal value of this parameter (e.g. we find 0(O) m ~ 0.92 for Uoo = 0.2) and that a second 
part of this branch of solutions exists, terminating at 0(0) = 0(O) m . It is labelled '2' on 
FigflQl When we decrease 0(0) on this part of the curve, we observe that the parameter 
b and the difference A = 0(oo) — 0(0) increases while value of parameter u\ decreases. 
However, the numerical results strongly indicate that A become very large (likely infinite) 
for 0(0) ~ 0.54 while b and u\ stay finite. The numerical analysis become very difficult in 
this limit. This suggests, however the occurrence of a critical value 0(O) C limiting the domain 
of existence of the second part of this branch in the parameter 0(0). For this choice of 0(0), 
the gap A becomes arbitrarily large and the dilaton function is very well approximated by 
0(r) = ln(r) in the intermediate region. 
One node for W 

The case of solutions with one node of the function W is also instructive and, we guess, 
leads to a pattern of solutions corresponding to several nodes of W. Increasing 0(0) from 
0(0) = — oo, we construct the first part of this branch of solutions corresponding to the 
deformation (by the field u(x)) of the first solution of the sequence of dilaton- Yang-Mills 
solutions constructed in [6| . This branch of one-node solutions can be constructed up to 
a maximal value of 0(0). For instance we find 0(O) m ~ —2.05 where we have ~ 0.126, 
b ~ 0.318, A w 2.5, as shown on FigfTTl This branch is labelled '1' on the figure. When 0(0) 
decreases on this part of the curve (which is equivalent to the limit — > 0), we observe 
that the minimum of W (denoted by Wm on the figure) progressively approaches —1 and 
the solution approaches the solution of [6] with W(0) = 1, W(oo) = — 1. 

However, this is not the end of the story; indeed the second part of this branch of solutions 
can be constructed for 0(0) < 0(O) m . This part of the curve is represented by the label '2' 
on FigJTTl Like the case of the zero-node solution, it seems that this curve is limited by a 
critical value of 0(0). Decreasing 0(0) on this part of one node branch (which is equivalent 
to increase in reference to Figs. [MID, we observe that the values of A and b increases and 
value of u\ decreases and that, again, A seems to go to infinity when 0(0) approaches the 
critical value 0(O) m . As a further check of this statement that the Bessel function plays a 
crucial role in the critical phenomenon, we report on Fig. [11] the numerical value of W m (the 
local minimum of W(r)) and observe that this quantity approaches the value W m ~ —0.16 
on the curve '2'. This coincides with the ratio W (x\) /W (x-i) given in Eq. (l29l . 
More nodes for W 

For solutions with n nodes of W, we expect that the scenario will be similar as the two 
cases we have analyzed in details. In particular, we conjecture that: (i) the second part of 
the curve corresponding to n— node branch of solution can be constructed up to a minimal 
value of the parameter 0(0); (ii) in the large Uoo-limit, the gauge function W of the non 
Abelian Euclidean solutions will approach the profile of the Bessel function between the 
n + 1-th extremum and infinity; (iii) in this same limit b will approach the absolute value 
of the second derivative of the function VFbessei at the corresponding extremum. Numerical 
investigation of the two-node solution supports this interpretation. A more involved check of 
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<j>(0) 

FIG. 10: The parameters A = (f>(oo) — (fi(0), b and u\ are plotted as functions of 4>(0) for = 
0.2 for zero node solutions. Curve A represents the values of the corresponding reduced action 
converted to the gauge </>(oo) = 0. 



this conjecture for generic n would require considerable numerical efforts with our numerical 
technique or a more refined analytical techniques of the type developed by D. Maison and 



collaborators 20, 21 



V. CONCLUDING REMARKS 



Motivated mainly from string theory we investigated the Euclidean version of a model 
involving a Yang-Mills field coupled to a dilaton field. We have constructed numerically 
different branches of finite action solutions of the resulting equations of motion. 

Euclidean theory is typically associated to tunnelling processes and to field theory at 
non-zero temperature. The solutions discussed in the present paper might be relevant to 
string theory at non-zero temperature. In any case they provide new saddle points in the 
Euclidean path integral. 

Our Euclidean model is very similar to static spherically symmetric sector of sponta- 
neously broken 577(2) gauge theory coupled to a dilaton field, which was investigated in 
22|, 23] and a discrete family of regular solutions were found corresponding to generalization 



of the 't Hooft-Polykov monopoles and their radial excitations. Main difference between Eu- 
clidean theory under investigation and static sector of spontaneously broken SU{2) gauge 
theory is in a way how dilaton field is coupled to electric part of the gauge field resp. Higgs 
field. 

The Yang-Mills-dilaton theory studied in the present paper can be viewed as a limit of 
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FIG. 11: The parameters A = <^(oo) — <^(0), b, u\ and the minimal value W m of the gauge function 
W are plotted as functions of ^(0) for = 0.2 for one node solutions. 



vanishing gravitational constant G — ► of the Einstein- Yang-Mills-dilaton theory investi- 
gated in 15[. Our results might help to find new branches of solutions in Einstein- Yang- 
Mills-dilaton theory. 

An interesting question concerning Euclidean solutions is the number of negative modes 
in the spectrum of small perturbations about these solutions. It is known that the instanton 
solutions, which describe mixing between equal energy states have no negative modes (at 
most have zero modes). Bounce solutions, which describe decay of a metastable vacuum have 



exactly one negative mode [24|, |25|, |26|, \2j\. Whereas presence of more then one negative 



mode typically shows that there are other solutions with lower action [28j]. For better 
understanding of the role of solutions which we discussed here and for their interpretation 
it is important to study question of negative modes about them. We plan to come to this 
topic in further investigation. 
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